Calculation of the current noise spectrum in mesoscopic transport: an efficient 

quantum master equation approach 

. JunYan Luo,^ Xin-Qi Li,^'Q and YiJing Yan^ 

' 'state Key Laboratory for Superlattices and Micro structures, Institute of Semiconductors, 

, ^ , Chinese Academy of Sciences, P.O. Box 912, Beijing 100083, China 

' ^Department of Chemistry, Hong Kong University of Science and Technology, Kowloon, Hong Kong 

(Dated: February 6, 2008) 

Based on our recent work on quantum transport [Li et al, Phys. Rev. B 71, 205304 (2005)], 
where the calculation of transport current by means of quantum master equation was presented, in 
Q ' this paper we show how an efficient calculation can be performed for the transport noise spectrum. 

, Compared to the longstanding classical rate equation or the recently proposed quantum trajectory 

method, the approach presented in this paper combines their respective advantages, i.e., it enables 
us to tackle both the many-body Coulomb interaction and quantum coherence on equal footing 
and under a wide range of setup circumstances. The practical performance and advantages are 
, illustrated by a number of examples, where besides the known results and new insights obtained in 

a transparent manner, we find that this alternative approach is much simpler than other well-known 
^ full quantum mechanical methods such as the Landauer-Biittiker scattering matrix theory and the 

pH , nonequilibrium Creen's function technique. 

o 

PACS numbers: 73.23.-b,73.63.-b,72.10.Bg,72.90.-f y 
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I. INTRODUCTION 

cn 
o 

\^ ' A modern trend in transoprt studies of mesoscopic systems is not only for the current-voltage characteristics, but also 
5 tor the noise properties QlSBBHHBaEa. The noise epect™, wh.ch .s a ~e ot the temporal co„el.t.o„ 

' . . . 

' between individual electron events, has been proved to be a unique tool to reveal different possible mechanisms which 
^ ' are not accessible by the mean current measurement In particular, via the measurement of current noise one is 
also able to extract information of the system parameters that govern the transport, as well as the internal energy 
Ch , scales of the mesoscopic system. 

^ , In principle, full quantum mechanical calculation for transport noise spectrum can be performed by either the 
Landauer-Biittiker scattering matrix approach [ill Il2l | or the nonequilibrium Green's function (nGF) technique [l^ 



Q. 



" [ijj. In practice, however, they were largely restricted to non- interacting systems. Very recently, the zero-frequency 
5^ noise of transport through quantum dot with strong Coulomb interaction was calculated [isf . with the nGF-based 
sophisticated real-time diagrammatic technique [l^. 

An alternative method to calculate transport noise is the classical rate equation approach 0,0,0] . This approach 
is much simpler than the full quantum mechanical methods, thus has been employed in some interesting transport 
systems, for instance, the important Coulomb blockade systems [iSI^- Nevertheless, owing to the classical nature of 
this approach, the quantum coherence that widely exists inside the transport systems cannot be described. Focused on 
the effect of quantum coherence in noise spectrum, Sun and Milburn studied the transport through a pair of coupled 
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quantum dots |23|. For this specific system, they derived a Lindblad-type master equation, then unravelled it to 
calculate the noise spectrum. With the same motivation to include quantum coherence , an alternative "n"-resolved 
quantum Bloch equation approach was proposed some time earlier by Gurvitz et al. |23l |. and recently applied it 
further (Here "n" stands for the electron number tunnelled through the junctions). We noticed that this 

approach started with the many-particle Schrodinger equation of the entire system, thus its applicability is restricted 
to zero temperature. Also, the derivation limited its validity condition in large bias voltage regime. In our recent 
work |25ll2a>l23j we extended this approach to finite temperature and arbitrary voltages, where we have based our 
derivation on the quantum master equation, which is free from any specific system, and can serve as a general and 
convenient starting point to study quantum transport. 

In Ref. 12 si we have established an explicit and compact expression for transport current and demonstrated its 
application by a number of typical examples. In this paper, we complete that master-equation transport formalism, 
by developing further the formulation for noise spectrum calculation. In next section, for completeness, we first outline 
the general idea of the "n" -resolved transport master equation, then present the compact formulas for transport current 
and noise spectrum. In Sec. IIIII we demonstrate the noise formula by the same examples in Ref. bsl i.e., transport 
through a non-interacting and interacting quantum dot, and a pair of coherently coupled quantum dots. For the 
simplest system of non-interacting dot, we easily recover the known result obtained first by Chen and Ting, which 
was based on a full quantum mechanical treatment in terms of non-equilibrium path- integral approach 13]. For 
interacting dot, we extend the result of zero- frequency noise obtained in Ref. to the whole frequency regime, and 
present brief discussion for the noise characteristics. For the coupled quantum dots, where the quantum coherence 
plays significant role, we show that our approach can be as efficient as the quantum trajectory method j^^. However, 
by noting that the quantum trajectory theory is based on the unravelling of Lindblad-type master equation, we believe 
that the formalism presented in this work is more powerful than the quantum trajectory approach, since not all the 
transport systems can be reduced to a Lindblad-type master equation (e.g. the above mentioned transport through 
interacting quantum dot). Finally, a brief conclusion is given in Sec. IIVI 



II. FORMALISM 



Consider the typical transport setup as schematically shown in Fig. 1, which consists of the central transport 
system (device) and two biased electrodes. The total system Hamiltonian reads 

H - Hsicl,Cf,)+ J2 Xl^"fe"afe""fe+ Yu Zl(*"'=f'^"fe^A' + H.c.). (1) 

a—L,B. k a—L,B. kfi 

Hs is the central system Hamiltonian, in which all the possible many-body interactions have been included, cj^ 
(c^i) is the electron creation (annihilation) operator of state "/i", which labels here both the orbital and spin states. 
The second term describes the non-interacting electrons of the two electrodes, and the third one describes timneling 
between the central system and the electrodes. 

For convenience, we re-express the tunneling Hamiltonian as 

^' = E(4^M + H.c), (2) 

where = Fl^ + Fr^ = jitakfj.o-afj.tT- By regarding the tunneling Hamiltonian as perturbation, the second- 

order cumulant expansion leads to the master equation for the reduced density matrix of the central transport system 

poo 

p{t) = -iCp{t) - dT{c'{t)g{t,T)c'{T)g\t,T))p{t). (3) 

Jo 
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FIG. 1: Schematic setup for electron transport through a mesoscopic system. 



Here the reduced density matrix p{t) is introduced by tracing over the electrode states from the density matrix px of 
the entire system, i.e., p{t) = l^rB[pT{i)\^ and (• • •) = TrB[(- • • )/9_b], with ps the density operator of the electrodes. 
The Liouvillian superoperators are defined as CA = [HstA\, C'A = [H',A], and Q{t,T)A = G{t,T)AG^{t,T), where 
G(t,T) is the usual propagator (Green's function) associated with only the system Hamiltonian Hs. 

The trace in Eq. Q is over all the electrode degrees of freedom, leading thus to the equation of motion of the 
unconditional reduced density matrix of the system. To describe the transport problem, we should keep track of the 
record of electron numbers that have tunnelled through the right (left) junction, in order to calculate the corresponding 
current through it. To this end, the average over states in the entire Hilbert space "B" in Eq. ^ should be replaced 
with states in the subspace which corresponds to n-electrons tunnelled through the right or left junction. Let 

us take the right junction as an example. The trace of Eq. (O over "_b("«)" leads to a master-type equation for the 
conditional reduced state p(""^(i) [2^ : 

-^tV^-'^^ct _ ct - 4->'"«-''ct - ct + H.c.}. (4) 

Here A^^^ = J2,yC!!^l,{±C)a^, and A^J^'' = J2a=L r^o^- The spectral functions C^t(±C) are defined in terms of 
the Fourier transform of the reservoir correlation functions, i.e., ci^^(±£) = J'^^dtc!^}j{t)e^^^* . The correlation 
functions in time domain are defined by {fi^{t)fav{T)) = ci^l{t - r), and {fati{t)fl^iT)) = Ci^llJ, - r), where (• • ■) 
stands for Ttb[{- ■ ■ )Pb''], having the usual meaning of thermal average. For p''"^\t), which describes the reduced 
system state conditioned by the electron numbers timnelled through the left junction, a similar equation as Eq. I^J can 
be obtained by interchanging the indices "L" and "_R". With the knowledge of the above conditional state p("^/^^(t), 
we are ready to calculate all the transport properties, such as the transport current and the noise spectrum. 



A. Transport Current 

Throughout this paper, we assume the unit system oi h = e = 1. Straightforwardly, the current through the ath 
(i.e. L or R) junction can be calculated via Ia{t) — d{Na{t)) / dt = X]n„ '^aTr[/6("")(f)], where Tr(- • • ) means trace 
over the internal states of the central system. Obviously, Tr[/9'"°^ (t)] = P{na,t) is nothing but the probability by 
that there are "no," electrons tunnelled through the ath junction until time "i" . Simple algebra leads to 

I^{t) = Tr[^j-VW+H.c.], (5a) 
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where the unconditional density matrix p{t), which is the sum of aU satisfies a simple unconditional master 

equation 

p{t) = -iCp(t) -[np{t) +Ti.c.], (6a) 
T^Pit) = IJ2 4"V(^) - Pmi+^] ■ (6b) 

B. Noise Spectrum 

Physically, the noise spectrum characterizes the temporal fluctuations of the entire circuit current. In steady 
state, the average currents through the left and right junctions are equal to each other, but the temporal (i.e. time- 
dependent) fluctuating currents are not. In general, the circuit current, which is typically the measured quantity in 
most experiments, is a superposition of the left and right currents, i.e., I{t) — alL{t) + bliAt). Here the coefficients a 
and b depend on the symmetry of the transport setup (e.g. the junction capacitances) and satisfy a + 5 = 1. In 
combination with the charge conservation law, say, II ~ Ir + Q, with Q the charge in the central system, we obtain 

I{t)I{Q) = alL{t)lL{Q) + blR{t)lR{Q) - abQit)QiO). (7) 

Accordingly, the noise spectrum is a sum of three parts 

S{uj) — aSL{oj) + bSR{uj) — abu^ Sq{uli)^ (8) 

where Sl/ji{ijj) is the noise spectrum of the current through the left (right) junction, and Sq{lu) characterizes the 
charge fluctuations in the central system. In the following we first develop general formulation for their calculation, 
then demonstrate a number of typical examples. 

For SL/fi{uj), we employ the MacDonald's formula [s^ 

Sa,{u})^2uoj dtsin{ujt) — {nl{t)), (9) 

where {nl{t)) = Y,^^ [Trp("°)(t)] = X;„„ nlPina,t). With the help of Eq. we can show that 

j^inlit)} = Tr[2^MiV"(t) + ^j+V + H.c.], (10) 

in which we introduce the "particle-number" matrix N°'{t) = naP^^^°'\t), and the superoperator means 

■ ■ ) = ^ E (• ■ ■ )4 ± 4(- ■ • )4;^] • (11) 
/J 

Then, S'c[(lj) is formally expressed as 

S^iiu) = 2c.Im[Tr{2(^j-)7V"(c.) + [^j-)7V"(-c.)]t) -t- (J;Wp(c.) + [^j+)p(-c.)]t) }] , (12) 

where N"{uj) = dtN°'{t)e^'^* , and p(cj) = dtp'^^e^^*' . Note that we are calculating the noise spectrum of steady 
state, thus the stationary density matrix p"* is used here. Simply, we have p{uj) = ip'^* ju). For 7V"(a;), we can first 
establish the equation of motion for iV" (t) , that is 

^Ar"(t) = -iCN'^if.) - [nN°'{t) - + H.c.]. (13) 

Then, from its Fourier transform 

i(w-/:)iv"(w) = 7^7V"(w)^-[7^iv"(-L^)]t-^j-)p(c^)-[^j-)p(-c^)]t, (m) 



5 



we can directly obtain N°'{uj). 

For the charge fluctuations on the central system, the symmetrized noise spectrum reads Sq{uj) = dT{N{T)N + 
NN{T))e"^^] = 4Re[/g°°dT5'(T)e''^^], where 5(r) = {N{t)N) = TvTvbP^ {t)NU{t)N p''^ pe], with U{t) = e-^"'' and 
N the electron-number operator of the central system. Using the cyclic property under trace, we have S{t) = 
Tr[iV (T(r)], where cr(r) = Ti b[U {t)N p'^'^ psU'^ {t)] is introduced. Noticeably, (T(r) is nothing but an alternative 
reduced matrix which satisfies the same equation of the usual reduced density matrix p(r), with the initial condition 
(t(0) = iVp"*. Straightforwardly, its Fourier transform a{uj) can be easily solved from 

i{LO - C)a{uj) = TZa{Lu) + [na{-Lu)]^ - N p'\ (15) 

Then, Sq{uj) — 'iY{,e{TT[Na{uo)\} . More specifically, in most cases we can carry out this quantity in the eigenstate 
basis of "TV" , that is 



^q(c^) = 4Re 



k 



(16) 



where Nk is the eigenvalue of N for the eigenstate and ak{i^) = {k\a(Lj)\k). 

To summarize, we have now constructed a general and quite compact formulation for quantum transport through 
mesoscopic systems, particularly for the calculation of noise spectrum which is usually a difhcult problem by any of 
other approaches. Obvious advantages of this formulation include its straightforwardness and simpleness, together 
with its applicability to handle the many-electron correlation, quantum coherence, and possible inelastic scattering 
(dissipative) processes during transport on equal footing. 



III. ILLUSTRATIVE APPLICATIONS 

In this section we demonstrate the above formulation by transport through first a non-interacting quantum dot, 
then interacting dot, and finally two coupled quantum dots. Remarks with respect to other approaches and completion 
of some results which are lacked in existing literature will be made in particular. 



A. Non-Interacting Quantum Dot 

For simplicity, we restrict to the simplest case that only a single dot-state is involved in transport process. In this 
example, the electron spin is an irrelevant degree of freedom which is thus neglected in the description (its inclusion 
only needs multiplying an entire factor 2 to the result obtained in the following). Accordingly, the dot Hamiltonian 
reads Hs = EqN , where N = c^c is the number operator. 

To proceed, we need to carry out the electron reservoir (i.e. the electrode) correlation function. Explicitly, 
cL^\t-T) = |t„PEfce=^*'^'=(*-"'Va^'(£^fe), and its Fourier transform ct\±C) = r„/i=^^(-£). Here, the wide-band 
approximation for the electrodes is applied, which results in the energy-independent coupling constant r„ — 2iTga\ta\'^ , 
where ga is the density of states, /i^'' is the usual Fermi function, and fa = 1 — fa- Since £"c = (— £'o)"c, we 
have A^a^ = clt\±C)c = Taf^\EQ)c. Neglecting the spin degree of freedom, only two states are involved, i.e., |0) 
and |1) for the empty the occupied states, respectively. With these identifications, the transport master equation © 
simply reads 



Po = -rLPo + Tijpi, 
Pi = -'^rPi+^lPo- 



(17a) 
(17b) 
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Hereafter, in order to carry out analytic expressions, we assume zero temperature. By inserting the stationary 
solution of Eq. H17|l into the current formula Eq. ISJ, the well known resonant current is obtained, i.e., Ir = —II = 
TlTrHVl + Tb). 

Now we turn to calculation of the noise spectrum. The current fluctuation spectra of the left and right electrodes 
follow Eq. (221l as 

Sl{^) = 2c^rLlm[-2iV(f(c^)+po(t^)], (18a) 
Sr{lu) = 2LoTRlni[2N^\Lu)+pi{Lo)\, (18b) 

where /5o('^) = ipf) /^i and pi(w) = ipf /ui, with pg* and pf the stationary solutions of master equation (|17|l . N^{uj) 
and N^{uj) can be evaluated based on Eq. ((T^ . For instance, Nq{uj) is given by 

zwiV(f(w) = rLN^{uj)-rRN^{cj), (19a) 
iiuN^iu;) = rRN[{ij)-TLN^{ij) + rLM^), (19b) 
and similarly, N^{lu) can be obtained. Then, straightforwardly, we have 



5l(w) = SRicj) = 21 



+ 



p2 p2 p2 



(20) 



where I = Ir — —Ij^ is the stationary current. Note that both the zero- and finite-frequency noise spectra of the two 
tunneling currents are identical for this simple system. 

To calculate the charge fluctuation spectrum on the central quantum dot, based on Eq. H16(l we simply have 
Sq{uj) — 4Re[CTi(a;)], where the fact that A^o = and Ni — 1 have been taken into account. Following Eq. H15(l . cti(w) 
is given by the following coupled equations 

iujao{uj) rLCTo(^) - rflCTi(cj), (21a) 
iiuaiiuj) = TRai{uj) -Tl(To{lu) - pi . (21b) 

Accordingly, we obtain Sq{u!) = 4//(r-^ -I- lu^). Combining all the three components of noise spectrum according to 
Eq. JSJ and assuming a symmetric conflguration (i.e. a = b = 1/2), we obtain 

aTlTr\ r2 



1+1 



(22) 



This is the well-known result for resonant tunneling through symmetric double-barrier structures, which was obtained 
first by Chen and Ting based on the non-equilibrium path- integral technique and also by Biittiker by using 
the scattering approach 0, where the so-called Fano factor F(w) = S{uj)/2I, especially its value at zero frequency 
F{0) = 1 — 2T]^Tr/T^, was discussed in particular. Here, we presented an alternative derivation which seems elegant 
and interesting. More importantly, its generalization to more complicated systems is straightforward, as we are going 
to illustrate soon. 



B. Interacting QD with Zeeman Splitting 



We now consider transport through an interacting quantum dot described by 

Hs^Y.^-^- + ^^^^i' (23) 

(7 

where Ey and are the electron levels with finite Zeeman splitting, and N„ — cJ.Co- is the electron number operator 
for spin a (it should not be confused with the reduced operator a(r) introduced in the calculation of SQ[iS) ). In this 
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model, the JJ-term accounts for the many-body Coulomb correlation. As in the previous example, we should first 
carry out the spectral function cj^\., , and the operator A^'' . Since the interaction does not flip the electron spin, the 
spectral function must be diagonal with respect to the spin indices, i.e., C^\, = S„„>Cc^a. Still under the wide-band 
approximation, we have Aq^'' = ^afa^\Wa)ca, where Wa- = E]'Sija + EiSia- + U{nijSia + 'T-i'^to-), which sets up four 
regimes for the applied bias voltage, i.e., (i) fiL > Ei > fiji, (ii) > E^,Ei > fiji, (iii) fiL> Ei + U,E^,Ei > fiR, 
and (iv) ^j,l> E^ + U, Ei + U, E^ , Ei> fiR. 

For the sake of being compact, here we only present detailed derivation for regime (ii). In this regime, the involved 
states include |0) (empty dot), 1 1 ) (occupation by a spin-up electron), and 1 1 ) (occupation by a spin-down electron). 
Accordingly, the transport master equation reads 



Pi = ~'^RPl + TlPo, 
Pi = -^F!.Pi + TlPo- 

Its stationary solution gives rise to the steady-state current I = Ir = —II = r/ {'2^ l + Ti?)- 
In the above mentioned representation of occupation states, SL/ji{Lj) are expressed as 

Sl{uj) = AuoTlI^[-2N^{uj)+p^{u)\, 

Sr{u) = 2cjrflIm[2iVf (cj) + 2iVf (cj) + p^{uj) + p|(cj)] . 



(24) 



(25a) 
(25b) 



Nq[uj), N^{u}) and N^{uj) can be obtained from Eq. (fTl|l . For Nq{ll;), for instance, the relevant equations explicitly 
read 



iu;N^{ij) - 2rLiVtf(t^)~rfl7Vf(c^)-rfl7Vf(a;), 
luN^^iuj) = TRNl'{Lu)~TLN^{Lu)+rLPoH, 
luNf^iuj) = TRNf:{Lu)-TLN^{Lu)+rLM^)- 



(26) 



Note that, in the above equations, p{uj) = ip'^* /uj, and p*** is the stationary solution of Eq. (|24(l . Straightforwardly, we 
obtain 



SLiio) = Sr{lu) - 21 



2i2TL)TR 



{2Tl + Tr)^ {2Tl + Tr^ {2Tl + Tr^ 



(27) 



Also, in the same occupation-state representation, Sq{lo) simply reads SqIlo) = 4Rc[(7^ {lu) +(T|(cli)], where cr-^ ^Lu) and 
cr|(a;) can be solved from 



iLuao{uj) = 2rLCTo(w) - r_R.(7|(cj) - r7jCT^(a;), 
iu!a^{uj) = r_R,fj|(w) - rLfJo(w) - pf, 
iujaiiu) = TRaiiu) ~ TL^oiuj) ~ pf. 



(28) 



Its solution gives 



AI 



{2Ti 



(29) 



In Table I and II all the components of the noise spectrum are listed, together with those of other voltage regimes, in 
which the derivation is completely the same as above. 

In Ref. the zero-frequency noise was carried out for the same system of this example, based on a complicated 
non-equilibrium real time diagrammatic technique. Here we extended the result to the whole frequency regime, using 
a much simpler and more straightforward approach. We notice that the noise spectra in regimes (i) and (iv) differ 
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TABLE L Zero-frequency and frequency-dependent noise components in the different bias regimes (i), (ii) and (iv), where 
r = r_L + r_R is the total level-broadening width. 



bias regime 


(i) 


(ii) 


(iv) 


S{0) 


2rLrfl(ri + r?;) 




4rirfl(ri-fr|) 




{2Vl + VrY 








ferirl 


srirl 




{2Vl+VrY (2ri, + rB)2+a;2 














{2Vl + VrY {2VL + rRY + Lo^ 


pa p2 _^ ^2 










r T^+Lo^ 


{2VL + rR) {2VL + rRY + Lu^ 


r r2 + (.j2 



TABLE IL Zero-frequency and frequency- dependent noise components in regimes (iii), where F = Fi, + Vr is the total level- 
broadening width. 



bias regime (iii) 



2VlTr{Vl + 2Vr) 3 2 s 

5(0) p5 -ii-Li-R] 

2rirfl(F3 - FiFfl + F^rj; + 5F|)a;^ 4Firfl(ri - 2r|)c^^ 

5'l(lj) F5(F2+t^2) + p3(p2_^^2)2 

6riFi;(Fz. + 2rfl) uj'' 4rir|(2F^ + 3Ffl) uo^ 

Sr{u}) p5 r^+uj^ (r2-i-tj2)2 

2 2FLrj;(r^ + 3F|) 4rirfl(FL + 2rfl) 

Sq p3 r2+(.j2+ p3 (F2+(j2)2 



from each other only by an overall factor 2. The reason is that the electron near the Fermi surface of the left electrode 
can only pass through the spin-down level of the quantum dot in regime (i), but can freely pass through the two levels 
in regime (iv) due to completely overcoming the Coulomb blockade. As a consequence, the result in regime (iv) is 
also the same as that of the non-interacting dot obtained in the previous subsection. In the Coulomb blockade regime 
(ii), electron can only pass through either the spin-down or the spin-up level in an exclusive manner. Compared 
to the situation of regime (i), the entrance probability of electron from the left electrode to the dot is enhanced by 
a factor of 2, whereas the leaving probability to the right electrode is the same. Therefore, in all the components 
of noise spectrum in regime (ii), "2rL" replaces the 'Tl" in the result of regime (i). In regime (iii), the Coulomb 
blockade is partially overcome, i.e., electron can enter the dot even there has been already an electron on the spin-up 
state in the quantum dot. The noise spectrum in this regime is relatively complicated, and seems beyond an intuitive 
simple interpretation. However, the unique features of noise spectrum in regime (ii) and (iii) may provide a pathway 
to distinguish Coulomb blockade phenomena from that resulting from non-interacting multi-levels, as emphasized in 
Ref.0 based on the zero-frequency noise. Moreover, the "asymmetry" of and Tji may provide a useful tool to 
determine their respective values, which is very important in molecular electronics, due to the need to identify the 
subtle connection of the molecule with the electrodes. 



C. Two Coupled Quantum Dots 



In this subsection we study the noise characteristics of transport through a pair of coupled quantum dots, as shown 
in Fig. 2. To highlight the underlying coherence effect, which surely goes beyond the scope of the classical rate 
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FIG. 2: Schematic illustration for transport through a pair of coupled quantum dots. 



equation approach, we would like to neglect the many-body Coulomb correlation. (Inclusion of it is straightforward 
but will make the solution more complicated). Accordingly, the system Hamiltonian reads 

Hs = Eic\ci + E2clc2 + n{c\c2 + c^ci). (30) 

To analytically carry out the operators A^^^ , it will be useful to diagonalize this Hamiltonian. To this end, we introduce 
a pair of new electronic operators, ci = uci + vc2 and C2 = uc2 — vci. The desired Hamiltonian Hs = £icjci + e2cjc2, 
gives rise to the diagonalization condition {E2 — Ei)uv + il(M^ + v^) — 0. Together with + ~ 1, we can obtain 
the values of u and u, and also the eigen-energies ei = Eiv? + i?2W^ + 2f7uu, and £2 — Eiv"^ + E2V? — 2Vluv. 

As done previously, under the wide-band approximation, we can first carry out the spectral function ci^\ then the 
operators A^^'' — T L[uf^j^\ei)ci — vf\^\e2)c2\, and = T ji[v f'"^\ei)ci + uf]j^\£2)c2\. Neglecting the irrelevant 
spin degree of freedom, the four relevant states of the double dots are |0) (empty dots), |1) (occupation of the left dot), 
|2) (occupation of the right dot), and \d) (occupation of both dots). In this state representation, Eq. jnj becomes 

Po = r_Rp2 - Tipo, 

Pi = Tipo + TflPd + «f^(/0l2 - /52l), 

h = -(Fl + r_R,)/92 - «f^(pi2 - /O21), (31) 

Pd = ^LP2 - ^RPd, 

P12 = iApri - ifl{pi ~ P2) - ^{Tl +^r)pi2, 

where p2i — Pi2^ ^^^d A = i?2 — -Ei- Note that the off-diagonal density matrix element pi2 purely arises from the 
quantum coherence between the two coupled dots, and has has no classical counterpart. From the stationary solution 
of Eq. (|31|l . we can calculate the steady-state current as 

J ^ ATlTr{Tl + Tr)^^ , . 

TLTj,[{VL + yRy + ^A^] + '^{yL + TR)^n^- ^ ' 



This result was also obtained in Ref. 

In the occupation-state representation, we express the electrode current fluctuation spectra as 



Sl{^^) = 2ujTLlni[-2Nl^[uj)-2N^[uj)+po{Lo) + P2{^)]. (33a) 
Sr{oj) = 2uoTRlm[2N^{w) + 2N;^{uj)+p2{uj) + pd{^% (33b) 
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where p{lu) = ip^* /uj, given by the stationary solution of Eq. (|31|l . For Sl{(^), Nq{u!) and N2{uj) are solved from 



(34) 



In the following we would like to present analytic result for A = 0, and numerical result for A 7^ 0. Also, for 
simplicity, we assume Tl = Tr = T. Straightforwardly, based on the solution of Eq. ||33J, we obtain S'l(ijj) = 
S'l(O) + + S'[{uj), where 



5l(0) = 



4m^(r^ -2r-^o^ + 8»-^) 

+ 4172)3 



(35a) 



2ri72 



r2 + 4f]2 r2 



(35b) 



2r3r]3(3r2 ~4»2) 

(r2 + 4f72)3 
r3r!2(i2r)2 _ r2) 



(r2 + 4f72)3 



T'^ + {uj- 217)2 
(w - 2U.)uj 



r2 + (t^ + 20)2 

{ijj + 2rj)tj 



(35c) 



_r2 + (w-20)2 r2 + (tJ + 20)2_ 

Note that 5'^(ct') is the well-known Lorentzian which stands for the incoherent component with similar structure as 
those shown in Table I and II, whereas S'l{uj) stems from the coherent coupling between the two quantum dots. For 
Sli{uj), precisely the same result can be carried out by the same procedures. For the charge fluctuations on the double 
dots, we have Sq{lu) — 4Re[ai{Lu) + ^2(0;) + 2cfdi^)], where 0^(0;) is obtained from 



iLL!ai{oj) 



-Tna2{uj) + rLCTo('^), 

-Tl'5^{uj) - Ti^ad{uj) - irj[CTi2(w) - CT2i(w)] - pf 

(Fl + Tu)'^2{io) + in[ai2{uj) - a2i(t^)] - P2*, 

-TLa2{^)+^R^d{Lo)-2p'^, 

-i^[ai(uj) - d2{oj)] + {Tl + rfl)5i2(cj)/2. 



(36) 



Then, we immediately arrive at 



Sq{uj) = 



8rn2 



(37) 



(r2 + 4f72)(r2 + t^2)- 

The noise characteristics is numerically shown in Fig. 3. In addition to Ref. some more discussions are briefly 
presented as follows: (i) Due to the coherent coupling between the two quantum dots, before the transport electron 
arrives at the right electrode, it may stay in the coupled dots for some time and oscillate between them. During this 
Rabi oscillation process, it prevents other electrons from entering the dots. As a result, a "dip" appears in the noise 
spectrum at the Rabi frequency. For example, for a setup parameterized by f2 — 4r, A = and a = b = 1/2, the dip 
locates at uj^i = 2fl, as shown by the solid curve in Fig. 3. This interpretation is illustrated further by the result from 
a pair of off-resonantly coupled quantum dots. Since for this case the oscillation frequency is cu'j^ = •\/A2 + (217)2, 
the "dip" should move to higher frequency. This is shown by the dashed curve in Fig. 3 where A = GF is assumed, 
(ii) In our approach, we explicitly distinguish the three components contributing to the noise spectrum, which enable 
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FIG. 3: Fano factor F{uj) versus the normalized frequency uj/T. The signature of quantum coherence between the two 
coupled quantum dots is the dip appeared at the Rabi frequency. In the numerical calculation, we assume Q = 4r. Other 
geometric parameters are labelled in the figure, and the corresponding results are shown by the solid, dashed, and dotted curves, 
respectively. Inset: zero frequency Fano factor, which shows a non-monotonic dependence on the scaled coupling strength il/F. 



us to gain deeper insight. For example, to demonstrate the effect of charge-number fluctuations on the coupled dots, 
we consider a very asymmetric configuration {a/b ^ 1). Noting that SL{i^) = [3, from Eq. © we conclude 

that the charge- number fluctuation noise Sq{uj) would play negligible role on the entire noise spectrum if a/6 S> 1. 
In particular, the difference of the noise spectra of a — b and a/b ^ 1 can give important information for Sq{uj). A 
simple comparison between the solid and dotted curves in Fig. 3 indicates that SQ{iLi) would reduce the entire noise 
spectrum. Also, for a/b ^ 1, at high frequency limit the noise spectrum approaches to unity, which is consistent 
with the result of transport through a quantum point contact [llj. (iii) Concerning the zero frequency noise [c.f. 
Eq. H35a|l ]. it shows a non-monotonic dependence on the (scaled) Rabit coupling strength fl/T. In the inset of Fig. 3, 
we plot the corresponding Fano factor F{0) = 5(0)/2/. Interestingly, it will be maximally suppressed at ~ 0.65r2/r. 
By increasing the coupling strength (51 ^ F), the zero- frequency Fano factor approaches 1/2, behaving similarly as 
the symmetric double-barrier resonant tunneling structure jllf . In the opposite limit, 11 <^ F, the Fano factor reaches 
unity which corresponds to the Schottky-type noise. 

Finally, we make a brief technical comment. As we have mentioned, this interesting system has been studied in 
Ref.H 

where the quantum trajectory approach was employed and particular attention was focused on the effect 
of quantum coherence. It was also proposed there that the quantum trajectory approach provides a powerful tool to 
calculate transport noise spectrum in the presence of internal quantum coherence in the transport systems. In this 
subsection, we have shown that the "n" -resolved master equation approach can easily solve those types of problems. 
Moreover, we notice that the quantum trajectory approach has severe limitation since it is based on the unravelling 
of Lindblad-type master equation. We believe that in many transport systems the simple Lindblad-type master 
equation cannot be obtained. A simple example is the one studied in the previous subsection. In this case, the 
quantum trajectory approach proposed in Ref. would fail. On the contrary, the "n" -resolved master equation 
approach can efficiently solve the problems in the presence of both quantum coherence and many-body Coulomb 
interactions. An interesting example may be the one studied in this subsection by including further the Coulomb 
interactions, which is currently under our study. 
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IV. CONCLUSION 

To summarize, based on the "n" -resolved quantum master equation we have presented an efficient approach to 
calculating of the current noise spectrum. It overcomes the drawbacks of the classical rate equation methods and the 
recently proposed quantum trajectory approach, since the former cannot account for quantum coherence, whereas 
the latter seems unable to handle many-body Coulomb correlations. In practice, this approach is much simpler 
and more straightforward than other well-known full quantum mechanical methods such as the Landaucr-Biittiker 
scattering matrix theory and the nonequilibrium Green's function technique. These advantages have been preliminarily 
illustrated by a number of examples of quantum transport through quantum dots, where the known results and new 
insights are obtained in a unified and transparent way. Further application to more complicated cases is worthwhile 
and straightforward, which will be the topic of our forthcoming work. 
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